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Abstract. White dwarfs in globular clusters offer additional possibilities to determine distances and ages of globular 
clusters, provided their spectral types and masses are known. We therefore started a project to obtain spectra 
of white dwarfs in the globular clusters NGC 6397 and NGC 6752. All observed white dwarfs show hydrogen- 
rich spectra and are therefore classified as DA. Analysing the multi-colour photometry of the white dwarfs in 
NGC 6752 yields an average gravity of logg = 7.84 and 0.53 Mq as the most probable average mass for globular 
CLi cluster white dwarfs. Using this average gravity we try to determine independent temperatures by fitting the white 

dwarf spectra. While the stellar parameters determined from spectroscopy and photometry usually agree within 
5_l ■ the mutual error bars, the low resolution and S/N of the spectra prevent us from setting constraints stronger than 

the ones derived from the photometry alone. For the same reasons the white dwarf spectra obtained for NGC 6397 
unfortunately do not provide an independent distance estimate of sufficient accuracy to distinguish between the 
long and short distance scale for globular clusters. 
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1. Introduction dwarf sequence as a standard candle for determining the 

distance to nearby globular clusters, similarly to the tra- 

As white dwarfs are the final stage in the evolution of j... , ■ c ,, , , , , , 

° ditional mam sequence fitting procedure using local subd- 
all low mass stars, many are expected to exist in glob- r i , • , • n T i-v.- u, 

, warts witn known trigonometric parallax. In tins case, trie 
ular clusters. However, due to their faintness and occur- , <. r , , , , 

' white dwarl sequence oi the cluster is compared to a se- 
rence in very crowded fields they were detected only af- , , , , , , c .,, , 

J J J qucncc constructed with local wliite dwarts with accurate 

ter the refurbishment of HST with the WFPC2. Several , . , . r™ +u a j-v, v a i. 

trigonometric parallax. Ihe method was then applied to 

candidate white dwarfs swere then soon identified inM15 NGC ^ (Renzini gt ^ M and 4? Tuc (Zoccali et 

(de Marchi & Paresce HSgBj , u> Cen (Elson et aLE©, al gjjjjj, where alg0 th regult from Renzini t &1 um 

NGC 6397 (Paresce et al. 1995a Cool et al. ESS, M4 , „„ co , , „ „ , JTT 

,„ , v , 1 — ■ _ ' , " . , for NGC 6752 was updated). However, while the updated 

fRichter et al. 1995, 2002). 47 Tuc (Paresce et al. 1995b -,. , t Hnr ,L„ „ .,, 

v , .,' — TfrL^ \^ . distance modulus tor NGC 6752 agrees well with all other 

Zoccah et al. EDQD andNGC 6752 (Renzini et al. ». digtance determinations for this c Lter, the white dwarf 

Renzini fc Fusi Pecci M suggested to use the white distance tQ 4? Tuc ig considerably shorter than that found 

Send offprint requests to: S. Moehler b y Gratton et al. HI997J lMi® via main sequence fitting. 

* Based on observations collected at the European Southern While it may seem strange to use the faintest objects in 

Observatory, Chile (ESO proposal 65.H-0531, 67.D-0201) a globular cluster to derive its distance white dwarfs offer 

** Based on observations with the NASA/ESA Hubble Space some advantages as standard candles when compared to 

Telescope obtained at the Space Telescope Science Institute, main sequence stars: 
which is operated by the Association of Universities for 

Research in Astronomy, Inc., under NASA contract NAS 5- - They come in just two varieties - either hydrogen- rich 
26555 (DA) or helium-rich (DB) - independent of their orig- 
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inal metallicity and, in both cases, their atmospheres 
are virtually free of metals. So, unlike in the case of 
main sequence fitting, there is not the problem of find- 
ing local calibrators with the same metallicity as the 
globular clusters. 

— White dwarfs are locally much more numerous than 
metal-poor main sequence stars and thus allow to de- 
fine a better reference sample. 

However, the method has its own specific problems, 
which are discussed in great detail in Zoccali et al. (2001) 
and Salaris et al. (2001). Indeed, the location of the white 
dwarf cooling sequence depends on: 

— the white dwarf mass 

On theoretical grounds, given the observed maximum 
luminosity reached on the asymptotic giant branch 
(AGB), the mass of currently forming white dwarfs 
in globular clusters should be 0.53 ± 0.02 M Q (Renzini 
& Fusi Pecci ll988l Renzini et al. H996|) . Unfortunately, 
there are no local white dwarfs in this mass range with 
directly determined masses (i.e. without using a mass- 
radius relationship). There is, however, a handful of 
local white dwarfs with spectroscopically determined 
masses near this value (cf. Table 1 in Zoccali et al. 
2001), which allows to construct a semi-empirical cool- 
ing sequence for Mwd=0.53M q , once relatively small 
mass-dependent corrections are applied to each local 
white dwarf. The spectroscopic determination of the 
mass of white dwarfs in a globular cluster was first 
attempted by Moehler et al. (|2000|l for white dwarfs 
in NGC 6397. However, the low S/N of the spectra of 
these very faint stars did not allow to determine the 
mass with sufficient accuracy. 

— the white dwarf envelope mass 

In the case of DA white dwarfs the cooling sequence 
location depends also on the mass of the residual 
hydrogen-rich envelope. This affects spectroscopically 
derived masses (see above), with the resulting mass 
being w 0.04 Mq higher when using the evolutionary 
envelope mass (~ 10~ 4 Mwd, Fontaine & Wesemael 
1997) as opposed to virtually zero envelope mass. This 
mass uncertainty corresponds to an uncertainty of 0™1 
in the distance modulus and 1-1.5 Gyr in the age de- 
rived from the main sequence turnoff. 

— spectral type 

DB stars are fainter than DA stars at a given colour, 
with the offset depending on the filter combination (i.e. 
the offset is greater in V vs. B — V than in / vs. V — 
I). However, more massive DA white dwarfs are also 
fainter at a given colour. 

The white dwarf sequence also allows to determine the 
age of a globular cluster if its faint end is detected. In that 
case one can derive the age of the globular cluster from the 
luminosity of its oldest and faintest white dwarfs. Aside 
from the observational difficulties and the uncertainties 
in the cooling tracks (see Chabrier et al. 120001 for more 
details) any error in the assumed mass affects the result. 



Table 1. Description of the used photometric data-sets. 
The filter name, the number of exposures and the total 
integration time (in seconds) are given. 



Filter 


No of exp 


Exp. time 


F336W 


12 


26,400 


F439W 


5 


10,000 


F555W 


5 


6,000 


F814W 


7 


7,000 



Recent very deep HST observations allowed to detect the 
white dwarf cooling sequence in M 4 to unprecedented 
depths of V«30 (Richer et al. I2002fl . As a preliminary re- 
sult Hansen et al. H2002[) derive an age of 12.7±0.7 Gyr 
from the white dwarf luminosity function of M 4, consis- 
tent with other independent age estimates (their error bar, 
however, does not include errors due to the uncertainty 
of the white dwarf mass) . Their result has recently been 
questioned by de Marchi et al. (20Q3|), who claim that the 
cluster membership of the white dwarfs can not be veri- 
fied down to sufficiently faint limits to obtain more than 
a lower limit of the age. Richer et al. ((2004 ) showed, how- 
ever, that different methods of data reduction and analy- 
sis account for the different depths reached with the same 
data and thereby defended the original result of Hansen 
et al. (UlU)- 

In view of the relevance of the white dwarf masses 
and spectral types to the problems described above we 
decided to observe spectra of white dwarfs in NGC 6397 
and NGC 6752 in order to determine their spectral types 
and get mass estimates. Pilot observations of the white 
dwarf candidates in NGC 6397 (63.H-0348) had already 
shown that the targets are hydrogen-rich DA white dwarfs 
(Moehler et al. 2000 ) , but did not allow much quantitative 
work. 

2. Observations and Data Reduction 

2.1. Photometry 

The dataset used for NGC 6752 in the present work con- 
sist of a series of HST/WFPC2 exposures taken through 
the filters F336W, F439W, F555W, F814W in March and 
April 1995 as a part of the HST programme G05439 (see 
Table P). Results based on this dataset have been already 
published in Renzini et al. I|1996)) and Ferraro et al. l|1997fl . 
The field is centered about 2' south-east of the cluster cen- 
ter (see Fig. 1 in Ferraro et al. I1997|l . 

All the photometric reductions have been carried out 
using ROMAFOT (Buonanno et al. HMty . a package 
specifically developed to perform accurate photometry in 
crowded fields. A version of ROMAFOT optimized to han- 
dle sub-sampled HST images has been used for the present 
work (see Buonanno & Iannicola 1988). Details of the re- 
duction procedure can be found in Ferraro et al. I|1997|l . In 
short, we combined all the images in each filter obtaining a 
median image free of cosmic ray events and other defects. 
In order to perform a deep search for faint blue objects 
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we used the combined F439W image (with an equivalent 
total exposure time of 10,000 sec) as the reference frame 
for object detection. The object list derived from the me- 
dian F439W image was used as input for the photometry 
of the median images in the other filters and preliminary 
colour-magnitude diagrams (CMDs) have been obtained. 
The preliminary CMDs have been used to select a list 
of objects that fell outside the main CMD sequences. A 
visual inspection revealed that most of them were spuri- 
ous objects (e.g., artifacts due to the diffraction patterns 
of bright saturated stars, hot pixels, etc). By using the 
cleaned object list, the standard PSF fitting procedure was 
finally performed on each single image by using a Moffat 
( 1969) function plus a numerical map of residuals in order 
to better account for the contribution of the stellar PSF 
wings. 

A final catalogue was then compiled with coordinates 
and instrumental magnitudes for each star in each filter. 
The magnitudes were calibrated to the HST flight system 
following Dolphin ( 2000). Specifically, we applied: 

1. The charge transfer efficiency (CTE) correction, which 
is negligible for the brightest white dwarfs, but reaches 
up to 0™25 for the faintest ones in the F336W filter 

2. The aperture corrections, computed by means of a 
sample of isolated stars in each WFPC2 chip and filter, 
in order to insure that the instrumental magnitudes re- 
fer to the counts in a / /5 radius around each star 

3. The gain correction, due to the fact that the four 
WFPC2 chips do not have identical response factors 

4. A correction to the U magnitude, removing the effect 
of the decreasing UV sensitivity of the CCD 

5. The flight system zero points, given in Table 6 of 
Dolphin PKJOjl 

The resulting colour-magnitude diagram is shown in 
Fig. H and clearly shows the white dwarf sequence. The 
photometric data for the twelve most isolated white dwarfs 
in NGC 6752, from which the spectroscopic targets were 
selected, are listed in Table |5] 

2.2. Spectroscopy 

The spectroscopic targets in NGC 6397 (65.H-0531(B)) 
were selected from the photometry by Cool et al. ©96), 
those in NGC 6752 (67.D-0201(B)) from the data dis- 
cussed in Sect. 12. ll We took care to select stars as isolated 
as possible in order to avoid contamination from neigh- 
bouring stars. The spectra were obtained in service mode 
with FORS1 at VLT-UT1 (Antu) and VLT-UT3 (Melipal) 
on the dates given in Table |3| We used the MOS mode of 
FORS1 (slit length 19") with the high resolution colli- 
mator (O'.'l /pixel), a slit width of ft' 8 and grism V300. 
The detector was a TK2048EB4-1 backside thinned CCD 
with 2048x2048 pixels of (24 ^m) 2 , which was read out 
with high gain (1.46 e~/count, 5.15 e~ read-out noise, 
NGC 6397) respectively low gain (2.87 e"/count, 6.07 e" 
read-out noise, NGC 6752) and normal read-out speed us- 
ing one read-out port only. This configuration yields a 
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Fig. 1. V vs U — V (Vega-based flight system) for the 
stars in NGC 6752. The stars listed in Tabled targets are 
marked by filled circles. 



reciprocal dispersion of 2.5 A/pixel and a wavelength cov- 
erage of at least 3500 A to 5700 A (depending on the 
position of the slitlet also larger coverage is possible) . 

The slit width is larger than the requested seeing of 
0'.'6 to allow for small positioning inaccuracies, but results 
in a seeing-dependent resolution in those cases where the 
actual seeing is smaller than the slit width (cf. Table[3J). As 
discussed below, due to the low resolution and S/N of our 
spectra the change in resolution due to varying seeing does 
not affect the results of the spectral fitting. As FORS1 is 
equipped with an atmospheric dispersion corrector MOS 
observations at higher airmass are not a problem. 

The spectra were supposed to be observed only if the 
seeing was better than 0'.'6 (corresponding to a resolution 
of 15 A). However, due to the faintness of our targets 
we had to use rather long exposure times (2955 sec for 
NGC 6397 and 4380 sec for NGC 6752), during which the 
seeing might have changed. A few exposures were stopped 
before reaching the requested exposure time due to tech- 
nical problems. The values provided in Tableware aver- 
aged over the exposure time unless noted otherwise. Some 
spectra were taken at too bad seeing (>1") or with rather 
bright and nearby moon. Those observations provided no 
useful data and are therefore not listed in Table |3| 

For NGC 6752 dome flat fields with two different illu- 
mination patterns were observed for each night. For the 
previous programme on NGC 6397, however, we had dome 
flat fields (again with two different illumination patterns) 
only for the nights of July 3, 8, and 24, 2000. For data 
from the other nights we used the flat field closest in 
time. As part of the standard calibration we were also pro- 
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Table 2. The photometric data for the twelve most isolated white dwarfs in NGC 6752 in the Vega-based HST flight 
system 



Star 


£^336 


*-»^336 


^439 




Vrrr 

►S55 


AVrrr 

* 555 


-^814 


*-* J 814 


WF9 941 


24™014 


U U4U 


9 c *™0 c i9 

tCi'J \JiJtCi 


o™090 


24™818 


o'TYno 

\J KjoyJ 


24™916 


U (JUL) 


WF2-3648 


22™919 


l;l 020 


24™152 


0™020 


24™119 


0™030 


24™195 


0™040 


WF2-5639 


21™707 


0™010 


23™261 


0™020 


23™281 


0™020 


23™447 


0™030 


WF3-296 


24™493 


0™080 


25™419 


0™020 


25™287 


0™040 


24™770 


0™070 


WF3-1610 


23™004 


0™030 


24™321 


0™020 


24™212 


0™020 


24™127 


0™030 


WF3-3584 


23™718 


0™040 


24™908 


0™020 


24™632 


0™020 


24™639 


0™040 


WF3-4909 


24™256 


0™070 


25™295 


0™020 


25™020 


0™040 


24™923 


0™050 


WF3-6441 


23™577 


0™040 


24™574 


0™020 


24" 1 480 


0™030 


24™613 


0™040 


WF3-6454 


23™750 


0™050 


24™895 


0™020 


24™628 


0™040 


24™266 


0™060 


WF4-1700 


22™983 


0™030 


24™210 


0™020 


24^100 


0™020 


24™067 


0™040 


WF4-2003 


23™733 


0™040 


24™730 


0™020 


24™450 


0™030 


24^577 


0™030 


WF4-2081 


23™397 


0^040 


24™520 


0™020 


24™370 


0™020 


24™522 


0™040 



Table 3. Observational parameters for the spectroscopic data. On August 10, 2001, the seeing monitor did not work 
and the seeing was estimated from the acquisition and through-slit images 



target 


start of observation 


exp. time 


seeing 


air mass 




moon 








[sec] 




distance 


illumination 


in vi v uoy 1 


9000/07/01 0^-1 8-94 01 1 


90^ 


u. oo 






007 


NGC 6397 


2000/07/01 04:09:25.965 


2955 


0"54 


1.148 


147? 4 


0.006 


NGC 6397 


2000/07/05 03:44:29.402 


2955 


0'.'72 


1.146 


115? 


0.153 


NGC 6397 


2000/07/05 04:35:30.895 


2955 


0'.'87 


1.174 


114? 5 


0.156 


NGC 6397 


2000/07/06 05:42:39.878 


2955 


0'.'80 


1.282 


102? 2 


0.254 


NGC 6397 


2000/07/06 06:33:39.752 


2955 


0'.'81 


1.425 


101?7 


0.259 


NGC 6397 


2000/07/08 03:49:31.331 


2955 


0756 


1.151 


80? 5 


0.453 1 


NGC 6397 


2000/07/08 04:40:30.861 


2955 


0"50 


1.192 


80? 1 


0.457 1 


NGC 6397 


2000/07/25 02:52:02.328 


2955 


0'.'53 


1.156 


124? 1 


0.422 2 


NGC 6397 


2000/07/25 03:43:02.144 


2955 


0'.'53 


1.204 


124? 5 


0.417 2 


NGC 6752 


2001/08/09 01:16:57.397 


1756 


07 67 


1.260 


86? 1 


0.796 3 


NGC 6752 


2001/08/09 01:53:59.110 


1724 


0761 


1.235 


86? 4 


0.794 3 


NGC 6752 


2001/08/10 00:57:05.415 


4380 


0763 


1.263 


95? 6 


0.715 4 


NGC 6752 


2001/08/10 02:16:43.582 


2400 


0"67 


1.229 


96? 2 


0.710 4 


NGC 6752 


2001/08/11 00:30:47.431 


4380 


0"65: 


1.287 


105? 


0.626 5 


NGC 6752 


2001/08/11 01:50:27.188 


4380 


0755: 


1.235 


105? 7 


0.620 5 


NGC 6752 


2001/08/11 03:18:27.263 


4380 


0750: 


1.280 


106? 3 


0.613 5 


NGC 6752 


2001/08/11 04:35:06.749 


4380 




1.418 


106? 9 


0.608 5 


NGC 6752 


2001/08/22 00:53:16.938 


4380 


0770 


1.237 


94? 6 


0.119 


NGC 6752 


2001/08/22 02:16:26.569 


4380 


0756 


1.262 


93? 9 


0.124 


NGC 6752 


2001/09/15 00:04:58.643 


4380 


0764 


1.239 


132? 8 


0.092 


NGC 6752 


2001/09/15 01:37:23.372 


4380 


0"60 


1.332 


132? 1 


0.086 


NGC 6752 


2001/09/18 00:25:56.746 


4380 


0765 


1.257 


100? 


0.006 


NGC 6752 


2001/09/18 01:41:59.797 


4380 


0"53 


1.363 


99? 3 


0.007 


NGC 6752 


2001/09/19 00:11:09.587 


4380 


0770 


1.250 


88? 2 


0.035 


1 moon rise 


at 11:58 UT 












2 moon rise 


at 00:47 UT 












3 moon rise 


at 02:05 UT 












4 moon rise 


at 02:56 UT 












5 moon rise 


at 03:49 UT 













vided with masterbias frames for our data. The masterbias 
showed no evidence for hot pixels and was smoothed with 
a 30x30 box filter to keep any possible large scale vari- 
ations, but erase noise. The flat fields were averaged for 
each night and bias-corrected by subtracting the smoothed 
masterbias of the night. From the flat fields we determined 
the limits of the slitlcts in spatial direction. Each slitlct 



was extracted and from there on treated like a long-slit 
spectrum. The flat fields were normalized along the dis- 
persion axis with a linear fit, which corrected most of the 
large scale variations. 

Due to the low resolution of the data we obtained a 
wavelength calibration spectrum with a slit width of 073 
(the center positions of the slitlets were unchanged) to 
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Fig. 2. Spatial brightness profile of the stars observed in 
NGC 6397 for each slitlet. The histogram profiles were ob- 
tained by averaging the two-dimensional spectra between 
4000 A and 5000 A. The straight lines mark the profile of 
the fitted sky background averaged over the same wave- 
length range. The triangles mark the spatial range over 
which the white dwarf spectra were extracted. 



allow a better definition of the line positions. The disper- 
sion relation for all spectra was obtained from this one 
calibration spectrum by fitting 3 rd or 4 th order polynomi- 
als along the dispersion axis. We used f 3 to f 6 unblended 
lines between 3600 A and 6200 A and achieved an r.m.s. 
error of about 0.1 A per CCD row. 

Due to the long exposure times the scientific observa- 
tions contained a large number of cosmic ray hits. Those 
were corrected with the algorithm described in Gossl & 
Riffeser piOlf . using a threshold of 15a and a FWHM of 
1.5 pixels for the cosmic ray hits. As the routine is not 
intended originally for the use with spectra we also re- 
duced the uncleaned frames to allow a check for any pos- 
sible artifacts of the cosmic cleaning procedure, but found 
none. The slitlets with the stellar spectra were extracted 
in the same way as the flat field and wavelength calibra- 
tion slitlets. The smoothed masterbias was subtracted and 
they were divided by the corresponding normalized flat 
fields, before they were rebinned to constant wavelength 
steps. 

For the NGC6397 data we then filtered the rebinned 
frames along the spatial axis with a median filter of 7 
pixels (corresponding to 0V7) width to erase any remain- 
ing cosmic ray hits and reduce the noise in the sky back- 
ground. We then identified regions uncontaminated by any 
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Fig. 3. Spatial brightness profile of the stars observed in 
NGC 6752 for each slitlet. The profiles were obtained by 
averaging the two-dimensional spectra (histograms) re- 
spectively the two-dimensional fitted background (smooth 
lines) between 4000 A and 5000 A. The slitlets containing 
WF3-6441, WF3-1610, and WF4-1700 clearly show the ef- 
fect of light from neighbouring stars providing a spatially 
varying background at the white dwarf's position. The tri- 
angles mark the spatial range over which the white dwarf 
spectra were extracted. 



stellar source in the filtered data (where the profiles of the 
stars are of course widened) and approximated the spatial 
distribution of the sky background by a constant or a lin- 
ear fit (cf. Fig. [21. The sky background spectra obtained 
this way from the filtered data were then subtracted from 
the unfiltered two-dimensional spectra. 

In the case of NGC 6752 the observed regions were 
in three cases (WF4-1700, WF3-1610, and WF3-6441) so 
crowded that the wings of spectra of brighter stars con- 
tributed to the background of the white dwarf spectra (cf. 
Fig. El ■ To account for this varying background we fitted 
the wings of the brighter stars with Lorentzian profiles. 
The FWHM and position of the two Lorentzian profiles 
were determined from averaged spatial profiles like those 
shown in Fig. [3] and the maximum was allowed to vary 
with wavelength to account for the spectral features of 
the bright stars. For WF4-2081 and WF4-2003 we per- 
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formed the sky background subtraction in the same way 
as described for NGC 6397. 



For both clusters the sky-subtracted spectra were ex- 
tracted using Home's (1986) algorithm as implemented in 
MIDAS (Munich Image Data Analysis System). Finally 
the spectra were corrected for atmospheric extinction us- 
ing the extinction coefficients for La Silla (Tiis I197?jl as 
implemented in MIDAS. As can be seen from Table EH we 
observed a large number of spectra for each target, which 
had to be co-added. Before averaging the spectra for each 
star we had to ensure that they are all on the same wave- 
length scale. As the wavelength calibration spectra are 
obtained during daytime slight shifts between the slitlet 
positions of scientific and calibration observations may oc- 
cur. Due to their very low S/N we cannot check for off- 
sets using the individual white dwarf spectra themselves. 
To account for possible zero-point shifts in wavelength we 
cross-correlated the range 5500 A-5650 A in the sky back- 
ground spectra of all scientific observations with the first 
spectra obtained on the night of June 30, 2000 (NGC 6397) 
and August 13, 2001 (NGC 6752). If the mean shift for all 
spectra from one exposure was more than 0.5 A (corre- 
sponding to 0.2 pixels) we corrected all spectra from that 
exposure with the mean offset. This correction however, 
cannot account for small offsets of the stars' positions from 
the center of the slit in dispersion direction. The signal-to- 
noise ratio of the individual spectra was much too low to 
obtain these offsets via radial velocities at a resolution of 
about 15 A. For the stars in NGC 6752 we were lucky to 
have a blue horizontal branch star in one of the slitlets. We 
used the spectrum of this star to determine any remain- 
ing velocity shifts and used those to correct the spectra of 
the white dwarfs to laboratory wavelengths. In NGC 6397 
we used the radial velocity offset derived from the spec- 
trum of WF4-358 for each observation to correct all other 
spectra. 



For a relative flux calibration we used response curves 
derived from spectra of LTT 6248, LTT 7987, and 
LTT 9239 with the data of Hamuy et al. lfT39"2]t for 
NGC 6752. For NGC 6397 we used the standard stars 
LTT 7379, Feige 110 (Hamuy et al. and CD 248 

(Oke 0.990). The response curves were fitted by splines 
and averaged for all nights. The white dwarf spectra of 
NGC 6752 were normalized with a curve fitted to the con- 
tinuum of the blue HB star spectrum, which had also been 
used for the determination of the velocity shifts. The re- 
sulting spectra are shown in Figs. [5] and H3 with the 
brightness of the stars decreasing from top to bottom. The 
region below 3800 A is completely dominated by noise, so 
that neither a spectral slope nor spectral lines can be iden- 
tified. This may be explained by the low UV throughput 
of FORS1, the reasons for which are not understood. Also 
the brightest star, WF4-358, however, shows no spectral 
lines below 3800 A. 




Fig. 4. Spectra of the three fainter white dwarfs 
NGC 6397 




'•COO 



lOOO 



4500 

x[A] 

Fig. 5. Spectrum of the brightest white dwarf observed in 
NGC 6397 



3. Fit of the photometry of NGC6752 

3.1. Theoretical magnitudes for the HST Photometry 

To calculate theoretical magnitudes in the so called 
"STMAG" system of the WFPC2 (Holtzman et al. ITMty 
we used the grid of DA model atmospheres of Koester 
(see Finley et al. 119971 for a detailed description) and the 
following definition of MgTMAG 



M s 



TMAG 



= -2.5 



J S(X)F X dX 



- 21.1 



JS(X)dX 

where S(X) is the filter response function and Fx is the 
mean intensity of the stellar disk or AHx in terms of the 
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4000 '1500 5000 



X[k] 

Fig. 6. Spectra of the white dwarfs in NGC 6752 (normal- 
ized) 

Eddington flux. Note that Mstmag is not an absolute 
magnitude but is used to distinguish theoretical from ob- 
served magnitudes to. With this definition of STMAG the 
magnitudes are all zero (Mstmag = 0) for a star with a 
constant F\ of 3.63 10" 9 ergs -1 cm -2 A -1 . 

A star observed in this system should have magnitude 
tistmag — 0, if the flux arriving at the earth is f\ = 
3.63 10~ 9 in the same units. The relation between observed 
and theoretical magnitude can be obtained from 

f ^ F 

ttR 2 

wstmag = -2.51og— J- + Mstmag + 21.1 

R 

= Mstmag + (m - M) - 5 log — + 41.991 

Here, R is the radius of the star, d the distance, and (to — 
M) the apparent distance modulus for the respective filter 
band. 

The observed magnitudes for the white dwarfs in the 
globular clusters were not determined in the STMAG sys- 
tem, but rather in the so called "WFPC2 Flight System" 
(Holtzman et al. H555l Dolphin I2UUU|) . We have chosen to 
transform the observations to the STMAG system by ap- 
plying the difference between the two zeropoints (STMAG 
from Holtzman et al. 115551 and FS from Dolphin 2000) 



before the fitting. The applied corrections (STMAG-FS) 
are 0™128, -0™724, -0™041, 1™228 for the magnitudes 
£^336 > -B439) V555, ^814 respectively. We have tested this pro- 
cedure in the following way: 

First we calculated the expected magnitudes for Vega 
in the STMAG system, using the absolutely calibrated 
Vega spectrophotometry from the HST archive (Colina et 
al. I1996fl . The inverse of the above corrections were then 
applied to the Vega STMAG magnitudes, which should 
bring them on the Flight System, with the resulting values 
-0^008, 0^076, 0^039, 0^002 for U 336l 5 439 , V 555l W 
The Flight System is defined in such a way that it should 
match the magnitudes of the ground based system used 
for pre-flight calibrations for stars with zero colors (see 
Holtzman et al. 1995 for a detailed discussion). Vega has 
approximately zero colors, but the magnitudes - though 
coming close - do not exactly match the observed UBVI 
magnitudes, which Holtzman et al. (1995) take as 0™02, 
0™02, 0™03, 0™035. One probable reason for this difference 
is that the calibration between flight system and ground 
system was built using stars (mainly white dwarfs) with a 
spectral type very different from that of Vega. The former, 
in fact, have much higher surface gravities and stronger 
Balmer lines. 

Although it would be possible to make additional cor- 
rections in an attempt to bring the system of the globular 
cluster white dwarfs closer to the Vega system, we do not 
believe that this is worthwhile in view of several other cal- 
ibration uncertainties and the observational errors for the 
faint objects. 

3.2. Fitting observations with theoretical magnitudes 

When trying to fit a set of observed magnitudes with the- 
oretical values we have in principle the possibility to de- 
termine simultaneously the three parameters T c ff, logg 
and distance. Effective temperature and surface gravity 
determine Mstmag as described above and provide the 
radius R of the white dwarf via the Wood |I995) ) mass- 
radius relation. In practice it turned out to be impossible 
to determine all three free parameters from only four ob- 
served magnitudes. In addition to observational errors, the 
major problem is the strong degeneracy with respect to 
the parameters. The relative energy distribution for the 
white dwarfs in the observed parameter range depends 
mostly on T e g and only very little on logg. The solid 
angle (irR 2 /d 2 ), which determines the difference between 
observed and theoretical magnitudes, depends strongly on 
both logg (through the mass-radius relation) and on the 
distance. Therefore a small change in logg can be easily 
compensated by a change in distance and vice versa. Fig.0 
shows the strong correlation between derived surface grav- 
ity and distance modulus. 

It was therefore necessary, to keep one or two param- 
eters fixed at some pre-determined value, to achieve a ro- 
bust fitting result. The fitting used a \ 2 routine based 
on the Levenberg-Marquardt algorithm (Press et al.D-992) 
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Fig. 7. Surface gravities determined from the photometric 
data listed in Table 0] for various distance moduli. The 
numbers refer to the position of the star in Table 2] 

and is described in more detail in Zuckerman et al. (2003). 
Before the fit, the magnitudes were dereddened assuming 
Eb-v = 0™04, and the extinction coefficients of Holtzman 
et al. 1(1335)1 . 

We decided to keep the distance modulus fixed as it 
is rather well known for NGC 6752: The distance modu- 
lus obtained by Renzini et al. (1996) using white dwarfs 
in the cluster was 13™05. A reevaluation of these white 
dwarf observations, using the "thick hydrogen envelope" 
models results in 13'"l5 (Zoccali et al. 2001). Recent deter- 
minations with main sequence fitting obtain 13™17 (Rcid 
and 13™23 (Gratton et al. 1TW7)> . Gratton et al. 
(2003) recently published a new distance determinations 
for NGC 6752 from main sequence fitting, which yielded 
(m-M)o = 13'"23 from Johnson photometry and (m-M)o = 
12-99 from Stromgren photometry. As it is unclear from 
their data which value is the correct one we do not use 
this new determination for our average distance. For our 
fits we have used the short distance originally determined 
by Renzini et al. I|1996l 13™05) as an extreme case, as 
well as the average of the more recent determinations not 
involving white dwarfs (13™20). 

The results listed in Table 0] show that the scatter 
in logg is quite small due to the very strong correla- 
tion between the radius (logg) and distance through the 
solid angle of the star (see above). The average value for 
the surface gravity is 7.96 for (m-M) = 13™05 and 7.84 
for (m-M)o = 13™20 with an error of the mean of 0.02. 
With the Wood (1995J) mass-radius relation for a typi- 
cal T eff of 15000 K and "thick hydrogen layer" (10~ 4 of 
the stellar mass) this corresponds to masses of 0.59 M Q 



and 0.53 Mq, for the short and long distance modulus, re- 
spectively. Assuming instead the more unlikely case of a 
"thin" hydrogen layer (< 10~ 6 Mwd) the mass would be 
0.56 M Q for (m-M)o = 13^05 and 0.50 M Q for (m-M) = 
13™20. If we consider that the progenitor stars in globular 
clusters should have less than 1 M Q and that NGC 6752 
has an exclusively blue horizontal branch with a very ex- 
tended blue tail, the smaller mass seems much more likely 
and gives strong support for the larger of the two dis- 
tances, in agreement with most distance determinations 
for NGC 6752. Using the distance modulus as variable 
parameter and logg as fixed (7.96, 7.84) yields essentially 
the same results, as expected. The effective temperature 
of the objects does not change significantly with the var- 
ious methods, and we may conclude that it is fairly well 
constrained. 

4. Spectral fits for NGC 6752 and NGC 6397 

The observed spectra for the white dwarfs in the two clus- 
ters were analysed with the same model atmosphere grid 
used for the calculation of the theoretical magnitudes. In 
order to verify the influence of different fit programs on 
the results (and thus get at least some estimate of the 
minimum true errors as opposed to the purely statistical 
errors provided by the fit routines) we used the proce- 
dure described in Finley et al. (1997, marked by "DK") 
and that described in Napiwotzki et al. ( 1999J marked by 
"RN"). 

The low S/N and resolution of the spectra leads to 
very large errors if T e ff and logg are both used as free 
parameters. We have therefore decided to keep log g fixed 
at the value 7.84 obtained as the most likely value from the 
photometry of NGC 6752. For the spectral fits we used a 
resolution of 15 A, corresponding to a seeing of 0'.'6 (which 
is the average seeing of our observations). 

To allow a better understanding of our results and their 
errors we discuss here some technical details of the fitting 
procedures as applied in this case. When searching for 
the best-fitting model (by the \ 2 fitting routine) in a first 
step the model is fitted to the observations at pre-selected 
"continuum" points. In the DK fitting routine these are 7 
points located in the continuum between H 7 and fig and 
longward of Ha in the red part of the spectrum, and in 
the middle between Balmer lines in the blue part. The 
observed continuum is determined from the median of the 
flux in a region 20 - 100 A wide around these points. 
For each of these continuum wavelengths a normalizing 
factor is calculated, which would adjust the model flux to 
the observations at these points. The complete normalized 
model spectrum is then obtained by interpolating these 
normalization factors. Since this uses a linear interpolation 
the normalized model can show abrupt changes in slope 
from one interval to the next if the observation is badly 
calibrated or extremely noisy, as is apparent at the blue 
ends of some of the spectra in Figs. |S]and|!|| After this 
normalization the x 2 is determined from the comparison of 
model and observations in several predetermined intervals, 
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Table 4. Fitting the photometry of the twelve most isolated white dwarfs in NGC6752 with T e ff and logg as free 
parameters, (m-M)o fixed at 13™05 and 13™20. Eb-v assumed as 0™04. 



Object 




(m-M)o 


= 13™05 






(m-M)o 


= 13™20 






T e a 


ATeff 


logg 


Alogg 


T c ff 


AT c ff 


logg 


Alogg 


WF2-241 


13000 


600 


8.08 


0.07 


13100 


500 


7.99 


0.06 


WF2-3648 


17400 


500 


7.94 


0.04 


17500 


400 


7.84 


0.04 


WF2-5639 


25200 


500 


7.86 


0.03 


25200 


500 


7.77 


0.03 


WF3-296 


9400 


1500 


7.66 


0.46 


9400 


1500 


7.54 


0.49 


WF3-1610 


16700 


1700 


7.94 


0.14 


16700 


1700 


7.84 


0.15 


WF3-3584 


13400 


1000 


8.00 


0.11 


13600 


1100 


7.91 


0.12 


WF3-4909 


11300 


200 


8.00 


0.02 


11200 


300 


7.90 


0.05 


WF3-6441 


14700 


1300 


7.99 


0.14 


14700 


1200 


7.89 


0.12 


WF3-6454 


10400 


1900 


7.52 


0.45 


10700 


2000 


7.47 


0.51 


WF4-1700 


16200 


900 


7.85 


0.08 


16300 


1000 


7.75 


0.09 


WF4-2003 


13000 


800 


7.86 


0.09 


13100 


700 


7.76 


0.09 


WF4-2081 


15300 


700 


7.98 


0.06 


15400 


600 


7.88 


0.05 


weighted mean 






7.96 


0.02 






7.84 


0.02 



which normally include the spectral line profiles. In the 
present case we have used the two intervals from 3870 - 
4480 and 4710 - 5100 A. These are the intervals plotted in 
the figures; the pure continuum region between 4480 and 
4710 A does not contain any useful information since the 
normalization always forces agreement. 

The major difference between the DK and the RN 
method is the definition of the continuum points. RN used 
two continuum regions of 40 A width adjacent to the fitted 
regions of the tFj, H 7 , and lines. As in the DK-case 
normalization factors between model and observed flux 
are computed and linearly interpolated for each individ- 
ual Balmer line. Since the high noise level combined with 
the strong overlapping of the lines makes the definition of 
continuum points between lines more and more difficult 
for the higher members of the Balmer series, the normal- 
ization of the region shortward of and including H e was 
done for all lines simultaneously. Four continuum regions 
(of 30 A width) were defined from the red wing of H e to 
3815 A. A linear fit of the resulting normalization factors 
was performed and applied to the fitted lines. 

Finally the best-fitting model is in both cases found by 
minimizing \ 2 a s a function of the model parameters log g 
and T c ff. This routine automatically determines errors for 
the fit parameters, which are however only statistical er- 
rors, assuming that the differences between the model and 
the observation are caused only by statistical measure- 
ment errors distributed normally around the correct value. 
There are, however, a number of sources of systematic er- 
rors, which cannot be determined easily. These include the 
reduction procedures, with errors in our case likely being 
dominated by the background subtraction. It is obvious 
that for an observation with a background both variable 
and several times higher than the continuum spectrum a 
small error in the background determination can have a 
profound influence on the line profile and thus the fitting 
results. Placing the continuum points differently, or inter- 
polating the factors quadratically instead of linearly will 
also change the results, but without prior knowledge of the 



exact observed spectrum there is no reason to prefer one 
method over the other. Such systematic errors can only be 
estimated by comparing results from different spectra for 
the same star, or comparing results from different authors 
and fitting/reduction methods. 

To verify the influence of the assumed resolution on 
the results we also performed fits for resolution of 10A 
and 20A, which yielded effective temperatures hotter re- 
spectively cooler by less than 1%, which is well below even 
the purely statistical errors given in Table|3] The observed 
spectra (confined to the ranges used for the model fits) and 
the best fitting models for NGC 6752 are shown in Fig.|SJ 

The effective temperatures derived for the white 
dwarfs in NGC6752 with different fit methods, although 
slightly different, usually agree within the errorbars (cf. 
Table and we therefore conclude that the choice of the 
fit methods does not affect our results. For WF3-6441 the 
spectroscopic temperatures agree with each other but dif- 
fer from the photometric one. This may be an effect of the 
rather problematic sky subtraction for this star (cf. Fig. [3] 
and Sect. I2.2|) . Despite our best efforts to approximate the 
spatial and spectral variation of the sky background in 
these crowded fields, we cannot exclude the possibility of 
residual background light in the extracted spectra. As the 
flux in the spectra is significantly lower than that of the 
background small errors in the background approximation 
can strongly affect the spectra. We therefore rather trust 
the photometric results for this star, despite its larger than 
average error. 

The observed spectra and best fitting models for the 
white dwarfs in NGC 6397 are shown in Fig. [5J For these 
stars we have only V and I photometry from Cool et al. 
( 1996J), from which we estimated temperatures in Moehler 
et al. H20001 cf. Table |SJ. These temperature estimates are 
very useful to distinguish between the hot and cool solu- 
tions obtained from the spectra alone (see TableHJJ). Based 
on the photometric temperature estimates we select the 
cool results of the spectroscopic fits for the three fainter 
stars and the hot result for the brightest star. As can be 
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Table 5. Fitting the white dwarf spectra in NGC 6752 with logg held fixed at 7.84. T e ff iP h is the result from the 
photometry for (m-M)o = 13*20. 



Object DK RN 





T c ff 


ATeff 


T e ff 


AT eff 


T e ff, p h 


ATeff^h 




[K] 


[K] 


[K] 


[K] 


[K] 


[K] 


WF3-1610 


15100 


900 


14400 1 


1000 


16700 


1700 


WF3-6441 


19200 


1000 


21600 


1300 


14700 


1200 


WF4-1700 


18500 


1600 


16800 


1800 


16300 


1000 


WF4-2003 


13400 


1700 


13400 


2700 


13100 


700 


WF4-2081 


12900 


800 


14600 


800 


15400 


600 


1 excluding 


the core 


of Hs 











Table 6. Fitting the white dwarf spectra in NGC 6397 with logg held fixed at 7.84. T e fjf )P h is the effective temperature 
estimated from VI photometry by Moehler et al. ( 2000 ) . We also give My from the model fit (assuming a thick hydrogen 
layer), the observed V corrected for a reddening of Eb-v = 0*18, and the derived distance modulus. 



Object 


V 


T e ff,ph 

[K] 


Teff 

[K] 


AT cff 

[K] 


DK 

My 


(m-M)o 


Teff 

[K] 


ATeff 

[K] 


RN 

My 


(m-M)o 


WF2-846 


23*76 


10300 


10700 


400 


11*73 


12 '"03 


11400 


700 


11*56 


12*20 








17900 


1200 


10*73 


13*03 


15700 


1000 


10*96 


12*80 


WF2-479 


23*33 


11500 


11500 


400 


11*54 


11 '"79 


12700 


800 


11*35 


11*98 








15300 


600 


11*01 


12*32 


14700 


800 


11*07 


12*26 


WF2-51 


23*46 


10800 


11200 


400 


11*60 


11*86 


12800 


1000 


11*34 


12*12 








15500 


500 


10*98 


12*48 


12800 


1000 


11*34 


12*12 


WF4-358 


22*19 


19400 


10200 


100 


11*88 


10*31 


10600 


100 


11*76 


10*43 








19000 


300 


10*62 


11*57 


17900 


300 


10*73 


11*46 



seen from Fig. the fitted model spectra for the hot (left 
panel) and cool (right panel) solutions are indistinguish- 
able by eye (except perhaps for the brightest star) and 
also the \ 2 values for both solutions are almost identical. 
As can be seen from Table |5] the cool solutions of the RN 
fit are hotter than those from the DK fit, whereas the sit- 
uation is just the opposite for the hot solutions, indicating 
the presence of systematic errors in our analysis. 

The parameters obtained from the spectral fit can be 
used to determine absolute magnitudes from the models, 
and - with the help of observed V values - also individual 
distances (see Table EJ. The average distance modulus for 
the three fainter objects is 11*9±0*1 [DK] and 12*1±0*1 
[RN], respectively (errors are statistical r.m.s errors only), 
yielding a mean distance modulus of 12*0, which is at the 
short end of the range of distance moduli found in the 
recent literature (11*99 - 12*24). 

The last object - although the brightest - is rather 
discrepant, and may be a foreground object or have a sig- 
nificantly smaller logg than the others. In view of the 
fact that the proper motions derived by King et al. I|1998f) 
support the cluster membership of all four white dwarfs 
analysed in NGC 6397 we decided to check the gravity 
for this star through a fit with both T e ff and log g as free 
parameters (see Table 0). Pilot spectra of this star with a 
combined exposure time of 1.5 hours yielded an effective 
temperature of 18,200±1300 K and a surface gravity of 
logg = 7.3±0.36 (Moehler et al. EHUDl- While the tem- 
peratures from all three fits agree very well, the scatter in 



surface gravity again shows that the systematic errors dis- 
cussed above are considerably larger than the statistical 
ones. The lower gravity obtained by the DK fit would sup- 
port the cluster membership of the brightest white dwarf, 
but would yield a mass of 0.48 M instead of the 0.53 M Q 
derived for the fainter white dwarfs in NGC 6752. 

5. Conclusions 

We observed a sample of white dwarfs in the globular clus- 
ters NGC 6397 and NGC 6752 and showed that they are 
all hydrogen-rich DA. From multicolour photometry we 
determined an average mass of 0.53±0.03M Q (uncertainty 
due to uncertainties of 0.02 dex in the average logg and of 
0*05 in the distance modulus of NGC 6752). This value 
(based on the assumption of a thick hydrogen layer) is 
identical to that assumed by Renzini et al. (1996) on the- 
oretical grounds. Therefore both observational and the- 
oretical arguments strongly advocate against the use of 
0.6 M Q (the mean mass of the local white dwarfs) or even 
more for the mass of hot white dwarfs in globular clus- 
ters when comparing observations to theoretical tracks. 
However, the limited S/N in combination with the low 
resolution of the spectroscopic data prevented the inde- 
pendent determination of masses from spectroscopic fits 
alone. Multi-colour photometry may be the better way to 
determine the physical parameters of white dwarfs in glob- 
ular clusters once their spectral types are known. For spec- 
troscopic observations our experience shows that crowding 
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Fig. 9. Spectral fits of the white dwarfs in NGC 6397. Only the part of the spectrum, which was actually fitted, is 
shown. The left panels show the hot solutions, the right panels show the cool solutions. 



can severely limit the usefulness of the data due to prob- 
lems with sky subtraction. We would therefore strongly 
recommend to look for white dwarf candidates in ground- 
based wide-field photometry to avoid the problems we en- 
countered for NGC 6752. Also a better sensitivity in the 
blue would be useful for future spectroscopy. 
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Table 7. Results for the brightest object in NGC 6397 
with logg as free parameter. We also give My from the 
model fit (assuming a thick hydrogen layer), the observed 
V corrected for a reddening of Eb-v = 0™18, and the 
derived distance modulus. 





T e ff 


AT cff 


logg 


A log g 


My 


(m-M)o 




[K] 


[K] 










DK 


18800 


340 


7.72 


0.06 


10" 1 47 


11™72 


RN 


18700 


520 


8.09 


0.12 


11™02 


11^17 
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Fig. 8. Spectral fits of the white dwarfs in NGC 6752. 
Only the part of the spectrum, which was actually fitted, 
is shown. 
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